Soliton ratchets induced by ac forces with harmonic mixing 
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The ratchet dynamics of a kink (topological soliton) of a dissipative sine-Gordon equation in the 
presence of ac forces with harmonic mixing (at least bi-harmonic) of zero mean is studied. The 
dependence of the kink mean velocity on system parameters is investigated numerically and the 
results are compared with a perturbation analysis based on a point particle representation of the 
soliton. We find that first order perturbative calculations lead to incomplete descriptions, due to the 
important role played by the soliton-phonon interaction in establishing the phenomenon. The role 
played by the temporal symmetry of the system in establishing soliton ratchets is also emphasized. 
In particular, we show the existence of an asymmetric internal mode on the kink profile which couples 
to the kink translational mode through the damping in the system. Effective soliton transport is 
achieved when the internal mode and the external force get phase locked. We find that for kinks 
driven by bi-harmonic drivers consisting of the superposition of a fundamental driver with its first 
odd harmonic, the transport arises only due to this internal mode mechanism, while for bi-harmonic 
drivers with even harmonic superposition, also a point-particle contribution to the drift velocity is 
present. The phenomenon is robust enough to survive the presence of thermal noise in the system 
and can lead to several interesting physical applications. 
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Transport phenomena based on nonlinear effects are at 
the heart of many problems in physics . In this context 
it was generally believed that an ac force of zero mean 
cannot lead to directed nonzero currents. Recent studies 
of the so called ratchet effect have shown that this belief 
was wrong A ratchet system can be described as a 
Brownian particle in an asymmetric periodic potential, 
moving in a specific direction in presence of damping, 
under the action of ac forces of zero average. The origin 
of a net motion is associated to the breaking of the space- 
temporal symmetries of the system leading to the 
de-symmetrization of Levy flights [for Hamiltonian sys- 
tems] 1^ and to phase locking phenomena between the 
particle motion and the external driving force [^j7| . 

This effect has a number of applications in various 
branches of physics and biology and is believed to be 
the basic mechanism for the functioning of biological mo- 
tors (see reviews Q and references therein). The ratchet 
effect, originally studied for Brownian particles, was gen- 
eralized to dynamical systems |^] and to partial differ- 
ential equations (PDE) of soliton type, mainly in the 
over-damped regime Q, or with asymmetric potentials 
l^^-Q. In the over-damped case, the damping in the sys- 
tem suppress all the degrees of freedom associated with 
the background radiation so that soliton ratchets become 
very similar to point particle ratchets. For under-damped 
or moderately damped systems, however, the situation 
is quite different since the radiation field can interact 
with the soliton and influence the transport. Under- 
damped soliton ratchets in asymmetric potentials and 
driven by sinusoidal forces, were recently investigated in 



Refs. [ |l3| , p^ . In particular, in Ref. the basic mech- 
anism underlying the phenomenon was identified in the 
existence of an asymmetric internal mode which couples, 
through the damping in the system, to the soliton trans- 
lational mode. Effective transport was found when the 
internal mode and the external force were phase locked. 
Moreover, it was shown that the effect of soliton trans- 
port decreases with increase of the damping, the maximal 
transport being achieved in the under-damped regime. 

On the other hand, it is known that for point par- 
ticle ratchets transport phenomena are possible also in 
symmetric periodic potentials, provided the driving force 
breaks suitable temporal symmetries of the system (i.e. 
the ones which relates orbits of opposite velocities in 
phase space fl^). Since in concrete applications it is 
more easy to act on the temporal part (by using external 
forces) than on the spatial part (by inducing distortion 
of the potential) of a system, it is interesting to explore 
the above possibility also in soliton systems. 

The present paper is just devoted to this problem. 
More precisely, we show that topological solitons of non- 
linear PDEs with symmetric potentials can acquire finite 
drift velocities in the presence of bi-harmonic forces of 
zero average consisting of the superposition of two har- 
monics, the fundamental and one of its overtones (har- 
monic mixing drivers). Bi-harmonic forces were used in 
the literature to suppress chaos in dynamical systems 
and in soliton equations as well as, to control the 
transport properties of single particle ratchets 0| . In the 
present paper we demonstrate, on the particular exam- 
ple of the sine-Gordon system, that bi-harmonic driv- 
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ing forces with certain symmetries can be effective to 
create soliton ratchets. The phenomenon is investigated 
both numerically, by direct simulations, and analytically, 
using symmetry considerations and soliton perturbation 
theory. We show that a first order perturbation analysis 
of the soliton dynamics captures only the qualitative fea- 
tures of the phenomenon, leading to poor quantitative 
agreements with numerical results. The reason of this 
discrepancy is ascribed to the soliton-phonon interaction 
which is obviously missing in a point-particle description 
(it arises only at the second order in the perturbation 
expansion ]17|). 

In general the situation can be described as follows. 
Besides a point particle contribution to the drift velocity 
there is an equally important contribution coming from 
the soliton-phonon interaction. This last manifests itself 
with the appearance of an internal oscillation on the soli- 
ton profile, asymmetric in space, which induces a net mo- 
tion in a similar manner as described in Ref. [p^ . In par- 
ticular we show that this oscillation can be phase locked 
to the external driving force and can couple to the kink 
translational mode, trough the damping in the system. 
The energy, "pumped", by the ac field into the inter- 
nal mode is then converted by the above mechanism into 
net motion of the kink. Internal oscillations on anti-kink 
profiles have opposite asymmetry compared to kinks, so 
that kink and anti-kink ratchets give rise to motion in 
opposite directions. 

The dependence of the phenomenon on system param- 
eters such as the damping in the system, frequency, am- 
plitudes, relative phase of the bi-harmonic driving force, 
as well as, on the presence of white noise in the system, 
is investigated by means of direct numerical integrations 
of the sine-Gordon equation. The interaction of soliton 
ratchets with the boundaries of a finite system is also in- 
vestigated. For reflecting edges we find that, depending 
on the initial velocity and position of the kink, the ratchet 
dynamics can be either destroyed or reflected as an anti- 
kink ratchet moving in the opposite direction. These 
results could be important for applications to physical 
systems such as long Josephson junctions, as we briefly 
discuss at the end of the paper. 

The paper is organized as follows. In Section I we 
investigate the dynamics of the perturbed sine-Gordon 
equation, as a model for soliton ratchets in presence of 
asymmetric forcing and damping, both in terms of sym- 
metry arguments and first order perturbation theory. In 
Section II we study soliton ratchets by direct numerical 
integrations of the perturbed sine-Gordon system and 
provide a consistent interpretation of the phenomenon. 
Qualitative and quantitative features of soliton ratchets 
are compared with the predictions of the first order per- 
turbation analysis. The phenomenon is investigated in 
the deterministic case (i.e. in absence of noise) and in 
presence of a white noise term in the system. We find 
that the soliton ratchets are robust enough to overcome 
the presence of small amplitude noises, this making them 
of interest for practical applications. In Section III the 



interaction of soliton ratchets with system boundaries is 
considered, while in Section IV we summarize the main 
result of the paper and discuss possible applications of 
the phenomenon. 

I. MODEL ANALYSIS 

Direct soliton motion induced by ac signals have been 
investigated in the literature mainly for ac forces with 
single-harmonic content . As well known, for symmet- 
ric field potentials this situation does not lead to soliton 
ratchet dynamics. To this regard we remark that the dc 
motion observed for a sine-Gordon kink driven by single 
harmonic forces in absence of damping |19|, and its gen- 
eralization to the case of small damping pO[, as well as, 
dc motion obtained from spatially-inhomogeneous drivers 
should not be confused with soliton ratchets. In 
these cases, indeed, the dc motion strongly depends on 
the initial conditions and quickly disappears as the damp- 
ing in the system is increased. On the contrary, soliton 
ratchets do not depend on initial conditions and exist 
also for relatively higher damping. We remark that net 
soliton motion independent on initial conditions, can be 
induced by the mixing of an additive and a parametric 
(single harmonic) driver as shown in Ref. ||2^ for solitons 
of the (f)'^ model. 

In this section we shall investigate soliton ratchets 
in symmetric potentials driven by periodic bi-harmonic 
forces of zero mean. As a working model we take the 
following perturbed sine-Gordon equation 

utt - Uxx + sinu =^ ~aut ~ E{t) + n{x,t) , (1) 

with a denoting the damping coefficient, n{x,t) a white 
noise term with autocorrelation 

< n{x, t)n{x', t') DS{x - x')6{t ~ t'), (2) 

and E{t) a driver of the form 

E{t) = El cos ujt + E2 cos {uiLjt + 9) , (3) 

(the case i?2 ^ is referred to as bi-harmonic driver with 
even or odd harmonic mixing, depending on m being even 
or odd). Note that although the symmetry properties of 
this driver are reduced if i?2 7^ and 9^0 mod tt, 
the force is periodic, with period T = 27r/w, and has 
zero mean (the analysis can be generalized to more har- 
monic components and to arbitrary non-sinusoidal peri- 
odic forces). The unperturbed version of Eq. (|^) [the 
perturbation being ef(t) = —E(t) — aut{x, t)] is the well 
known sine-Gordon equation with exact soliton (kinks, 
anti-kinks) solutions which depend on a free parame- 
ter, the velocity v of the kink, which lies in the range 
-1 < u < 1. 

It is of interest to investigate the conditions under 
which a bi-harmonic driver of type (^) can induce soli- 
ton ratchets in Eq. (Q). To this end we remark that due 
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to the translational invariance of the sine-Gordon system 
(we assume an infinite system or a finite one with peri- 
odic boundary conditions) we have that to each soHton 
trajectory with velocity v there is a specular trajectory 
with the velocity —v. One can expect that, in analogy 
with single particle ratchets jl^ , only forces which break 
the V — *■ —V symmetry should induce net motion. This 
argument can be formalized in terms of the kink velocity 
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xUxtdx 



(4) 



as follows (note that one could use the momentum of 
the kink as well, instead of the velocity). Among the 
possible shifts and reflections in t, x and w, we identify the 
symmetry operations which change the sign of v keeping 
the sign of the topological charge 



Q = 



2tt 



Uxdx , 



(5) 



unchanged (this means that we avoid kink anti-kink 
transformations). It is easy to check that there is only one 
symmetry transformation which changes the sign of kink 
velocity and leaves the equation of motion unchanged, 
i.e. 



X —X -l-xo, t ^ t + — , u ^ —u + 27r . (6) 

This holds true for driving fields satisfying the following 
condition 



E{t + r/2) = ^E{t) 



(7) 



We remark that the above symmetry argument accounts 
only for the breakage of the v — > —v point-particle sym- 
metry of the soliton, ignoring possible contributions com- 
ing from the soliton- phonon interaction p^ . From Eq. 
(^ we see that condition (0) is always satisfied for drivers 
with odd mixing while it is always broken for drivers with 
m even (obviously we take E2 ^ Q). Thus, the v — s- —v 
symmetry predicts that a sine- Gordon soliton should ex- 
hibit a ratchet dynamics when driven by a m = 2 force 
(one needs to break the symmetry (|^) to get the drift 
motion), but not when driven by a m = 3 force. 

To the same conclusion one can arrive also from first 
order perturbation theory, taking as collective coordinate 
ansatz for the kink profile 



w(x, t) = 4 arctan exp 



X{t) 



(8) 



An ordinary differential equation (ODE) for kink's center 
of mass X{t), can be readily obtained by differentiating 
the momentum (see Refs. |l9| , p^ ). 



-t-00 



fdx 



(9) 



with respect to time, and using Eq. (|l|) and ansatz (||) 
to simplify the expression. This leads to 



dP_ 
'dt 



= -aP + 2TTE{t) , 



or equivalently 
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■i) = --(l-w^) -TTy/l - v'^Eit) + 4a?; 



(10) 



(11) 



where v{t) — X{t) and we assumed the usual relativistic 
relation P{t) = 8w(t)/-\/l — v{tY between velocity and 
momentum to be valid for all times. Equation (HOh can 
be readily solved for P(i), from which the kink velocity 
can be obtained as 



P{t) 



8a/1 



P2(t) 



(12) 



From this equation an analytical expression of the aver- 
age kink velocity 



— / v{t')dt' , T = — , 



(13) 



valid in the limit Ej/ ^ o? + (j<^)^ << I7 3 = li2, (i.e. 
small momentums or small drift velocities), can be ob- 
tained for the case of m = 2, by expanding the square 
root in Eq. (n2) in series, this giving 



EIE^-k'^ 



512 (a 



arctan 



2 -t-a;2)Vc? + 4cj2 
aip? 3^2) 



i(e - ^o) 



(14) 
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From this expression we see that the dependence of (w) 
on the relative phase is perfectly sinusoidal. Similar cal- 
culations for the m = 3 case show that the average kink 
velocity is zero independently on the value of the relative 
phase Q (as well as of an arbitrary initial phase). 

These results can be easily understood from the effects 
of the symmetry properties of the force on the dynamics. 
In Fig. 0the force -B(t), viewed as sequence of alternating 
pulses of equal intensities (i.e. with the same area under 
the curve) and indicated by dark and light fillings in the 
figure, is reported for the cases m = 2,3. We see that 
for m = 3 these pulses perfectly balance so that no net 
motion can arise, while for m = 2 there is not such a 
balance (in both cases, however, the average of the force 
is zero). 

In Fig. ^ we report the dynamics obtained from nu- 
merical integrations of Eq. ( pT| ) for the case m = 2. We 
also show the response of the system to the single dark 
and light pulses composing the force, from which we see 
that although these pulses have equal intensities, the an- 
swer of the system is quite different in the two cases. 
Note that the negative pulse is more effective than the 
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positive one to produce motion as one can see from the 
fact that the area under the negative trajectory is greater 
than the one generated by the positive pulse. This is ob- 
viously a consequence of the nonlinearity of the system 
(in a linear system the area under the two curves would 
just be the same). For the particular example of Fig. ^ 
(i.e. 9 — 0) one expects then that a net motion in the 
negative direction will exist . This is indeed what one 
finds from integrations of the perturbation equation ( pT| ) 
with e = [see Fig. (|) below]. 

In Fig. a the motion of the kink center of mass, X , 
obtained from numerical simulations of Eq. (11) for the 
cases m = 2, 3, is also reported. We see that while a 
well defined drift velocity in the m — 2 case arises, no 
dc motion is present in the m — S case, this being in 
agreement with our symmetry analysis. 

We also find that the average velocity in (^^, com- 
puted after the system reached steady state regime, de- 
pends on the relative phase 9 with a low which is be well 
approximated by (v) — Asm{d — Oq), with A = 0.058 and 
6*0 = 0.8. This is shown in Fig. ^ where the dependence 
of the average velocity on 9, as computed from Eq.(pT|), is 
reported with dots, while the continuous line represents 
the above approximating function. Note that, although 
the sinusoidal dependence is in perfect agreement with 
the approximate result in Eq. (^, the explicit values 
of A, 9o differ from those predicted by ( [l^ ) , these being 
respectively A = 0.12, 9o = 0.56 (this is due to the fact 
that, for the chosen set of parameters, the approximation 
of small drift velocities is not valid). 

We also note that from Eq. ( p^ ) the point-particle con- 
tribution to the drift velocity is expected to be cubic in 
the driver amplitudes, this implying that there could be 
equally important contributions to the soliton velocity at 
higher orders of the perturbation expansion. In partic- 
ular, the interaction of the soliton with the phonons in 
the system, first appearing at second order ||lj], should 
not be overlooked. Although the development of a theory 
which includes second order effects is quite challenging, 
it is out of the purposes of the present paper. In the next 
section will shall instead resort to numerical simulations 
of Eq. (P for a full investigation of the problem and 
compare the results with those predicted by the present 
section. 



II. NUMERICAL STUDIES 

To numerically investigate sine-Gordon soliton ratch- 
ets driven by bi-harmonic fields, we have used standard 
finite difference schemes to reduce Eq. (^ into a set of 
ODE which were then integrated in time with a 4th order 
Runge-Kutta method. To understand the basic mecha- 
nism underlying the phenomenon we shall first concen- 
trate on the deterministic case by putting n{x,t) = 
in Eq. (|l]), and then show that the results obtained in 
this case will survive in the presence of noise. In Fig. 



H the dynamics of a sine-Gordon kink, initially at rest, 
driven by a bi-harmonic driver with m = 2 and by a 
single harmonic driver (i.e. E2 — 0), are reported in 
Figs. Ha,b, respectively (note that we use contour plots 
to show the time evolution surface generated by the kink 
profile). From these figures we see that, in the case m — 2 
the soliton center of mass move with a constant drift ve- 
locity (note that the shape of the kink during the mo- 
tion is highly distorted), while for the single harmonic 
driver, it oscillates around the initial position (periodic 
boundary conditions are used in the simulation). This 
demonstrates the importance of bi-harmonic drivers in 
establishing soliton ratchets. 

In Fig. ^ we report the dynamics of an anti-kink 
ratchet, for the same parameters values as in Fig. |^. 
We see that, in analogy with soliton ratchets in spatially 
asymmetric potentials [^,|l0|,^ , anti- kinks ratchets move 
opposite to kinks, the absolute value of the drift velocity 
being the same. 

To check the point particle perturbation analysis of 
the previous section, we have studied the dependence of 
kink's average velocity on system parameters for the case 
of a bi-harmonic driver with m — 2. 

In Fig. ^ we report (w) as a function of a for differ- 
ent values of system parameters. Note that the curves 
display similar behaviors, with a resonance peak in the 
under-damped regime and a quick decay to zero at higher 
damping. This behavior is similar to the one reported in 
Ref. ||l^ for soliton ratchets in asymmetric potentials, 
and suggests a possible common mechanism of the phe- 
nomenon (see below) . We also remark that the interrup- 
tions of the curves (cut-offs) at small a values are due 
to the disappearance of the kink as a consequence of the 
onset of spatio-temporal chaos in the system (the cut- 
offs delimit the borders of the existence diagram of the 
kinks). 

Panel (a) of Fig. ^ shows the dependence of {v) on a 
for different values of the driving frequency. Note that 
the velocity is influenced by resonances with the plasma 
frequency and its harmonics, as one can see from the non- 
monotonous behavior of the cut-offs at small a (the cut- 
off velocities in this case have their largest values close 
to a; = 0.5 and to = 1). 

Another interesting property emerging from this figure 
is that the kink velocity is enhanced at low frequencies, 
and the average velocity decreases by increasing the fre- 
quency (already for frequencies above u = \ directed 
motion is hardly visible). This is a consequence of the 
kink inertia to react to fast oscillations. On the contrary 
drift velocities are observable also at quite small values of 
uj (at very small values, however, the dynamics becomes 
complicated and requires long computational times - we 
checked explicitly the case uj = 0.01 for which an aver- 
age drift is still visible). The dashed lines in the figure 
represent the results of the point particle perturbation 
analysis of the previous section. We see that the agree- 
ment, although qualitatively reasonable, is not so good 
from a quantitative point of view. 
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In panel (b) of Fig. |^ the dependence of (v) on a is 
reported for different driver amplitudes (for simplicity we 
have varied Ei, fixing the ratio E2/E1 — 0.65 ). We ob- 
serve a situation similar to the one shown in Fig. |^a. 
By increasing the driver amplitude the system reaches 
the chaotic regime and cut-off values in a quickly ap- 
pear. The resonant peak in this case is quite weak and 
visible only for some narrow range of driving amplitudes. 
Also here the predictions of the point particle perturba- 
tion theory are quantitatively quite poor. In the panel 
(c) of Fig. 7| we have shown the dependence (v) (a) for 
different va ues of the relative phase 0. We see that by 
changing 9 one can change a maximum of the curve at a 
given value of a, into a minimum. This is a consequence 
of the sinusoidal dependence of the average velocity on 
6 predicted in the previous section. To show this, we 
have reported in the inset of the figure (v) versus 9 for 
a fixed value of a, as computed from direct numerical 
integrations of the sine-Gordon system. We see that the 
numerical points are well fitted by a sinusoidal law as 
expected from the first order perturbation result of the 
previous section. We remark that a similar sinusoidal 
dependence was also found in Ref. [Q for single ODE 
ratchets, this confirming the existence of a point-particle 
contribution to the effect. In experimental situations in 
which the relative phase between the two drivers is not 
accessible, one should consider 9 a,s a. random variable, 
and a final average on it should be taken. The above re- 
sults then imply that no drift velocity can exist in these 
cases (soliton ratchets can be induced only if the relative 
phase 9 remains constant in time). 

It is also interesting to note from Fig. ^ that reversal 
currents can be induced by changing the relative phase. 
In Fig. ^ we show how the curve 9 = tt of Fig. (which 
display current reversal at low damping), changes as the 
driving amplitude is increased. We see that by increasing 
the amplitude of the driver the kink velocity is increased, 
this leading to an upwards shift of the curve. This means 
that current reversal observed for some value of 9 can 
be removed by properly adjusting the driving amplitude 
(and vice- versa) . Note that a further increase of the am- 
plitude can change the shape of the curve destroying the 
resonant-like character. 

In Fig. ^ we report the dependence of the average ve- 
locity on the driver amplitude Ei with the ratio E2/E1 
fixed to 0.65, and with the relative phase between the 
two drivers 9 = 1.61. From this figure it is clear that 
the kink drift velocity depends non linearly on the driver 
amplitude Ei, with an almost cubic low as one can see 
from the log — log plot in the inset. This result is in 
good qualitative agreement with our perturbation anal- 
ysis (note that the EfE2 dependence in Eq. (^4|) implies 
a Ef dependence if one fixes the ratio E1/E2, as done 
in the numerical simulations). We remark, however, that 
the case uj = O.II, denoted by * in Fig. |[ indicates a 
deviation from this law at higher values of Ei . 

We also checked the predictions of point particle sym- 
metry arguments and perturbation theory for the case of 



bi-harmonic forces with odd harmonic mixing. In Fig. 
[lo| the dynamics of a kink driven by a bi-harmonic force 
with m = 3 is reported. In contrast with the prediction 
of the previous section, we see that kink can acquire a 
drift velocity also in this case. The direction of the mo- 
tion, however, depends not only on the relative phase, 
but also on the initial phase (or initial time) of the driv- 
ing force. By changing the initial phase we can achieve 
soliton moving in the opposite direction with the same 
velocity, thus by averaging on the initial phase one gets 
a zero mean velocity for the kink [simulations of the kink 
dynamics with different initial phases in the interval [0, T[ 
show that roughly half of that interval of initial points 
yield attraction to kink solutions moving to the right, 
the other half leading to kink solutions moving to the 
left with the same velocity] . A similar study for the case 
of m = 2 did not show any dependence on the initial 
phase or on initial time. Due to the sensitivity on initial 
conditions, we conclude that the kink motion in the case 
TO = 3 is not related to the ratchet phenomenon (this is 
similar to the cases reported in Refs. I9,2C|]). 

The net motion observed for the m = 3 case for fixed 
values of the initial phase, however, is an interesting phe- 
nomenon to explore by itself, since it is not linked with 
point-particle features of the soliton dynamics (these are 
excluded by the results of the first order perturbation 
theory), and is entirely related to the soliton-phonon in- 
teraction in a similar way as discussed in Ref. |p^ . 

From the above analysis the following conclusions can 
be drawn. Although first order perturbation theory cap- 
tures some qualitative feature of soliton ratchets, it does 
not provide a satisfactory description of the phenomenon. 
This is clear both from the fact that there is a poor quan- 
titative agreement between the PDE results and first or- 
der perturbation analysis in the case of to = 2, and from 
the fact that for to = 3 it fails to predict the existence 
of a drift velocity depending on an initial phase. The 
reason for this discrepancy is that this analysis includes 
only point-particle aspects of the soliton dynamics, ig- 
noring completely the internal structure of the soliton. 
In analogy with the mechanism described in Ref. [p^ , 
one could expect a strong contribution to transport com- 
ing from the soliton-phonon interaction. Since this in- 
teraction arises only at second order in a perturbation 
expansion, this explains why first order calculations fail 
to capture the phenomenon. 

To elucidate the mechanism underlying soliton ratch- 
ets it is useful to investigate the kink dc motion in more 
details. In Fig. ^ we depict the kink profile while exe- 
cuting the ratchet dynamics in the case of a m = 2 driver. 
From this figure the existence of an internal oscillation 
(internal mode) of the kink profile is clearly seen. Note 
that the internal mode is strongly de-symmetrized with 
respect to the center of the kink . In the panel (a) the 
kink moves from the left to the right and the mode ap- 
pears behind the kink. We checked that the motion of 
the kink is locked to the external drive. This can be seen 
from Fig. pT|b, which shows the kink dynamics during 
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two periods of the driving force E{t). This is in agree- 
ment with the results of Ref. The last panel, (c), 
shows the profile of the kink while moving in opposite 
direction than the one in panel (a). Note that when the 
kink moves to the left, the internal mode is on the right 
side from the kink center, so it is again behind the kink. 
This result indicates that there is an evident contribu- 
tion to the directed kink motion hidden in the asymmet- 
ric character of the internal mode and in its interaction 
with the kink center of mass. We also checked that by in- 
creasing the damping the internal mode become smaller 
and smaller and almost disappears in the over-damped 
limit. This correlates with the numerical results showing 
the mean velocity of the kink rapidly decreasing with the 
increasing of the damping. 

It is worth to note that although the are no internal 
modes frequencies in the spectrum of the small oscilla- 
tion problem around exact soliton solutions of the pure 
sine-Gordon equation, these can arise from the pertur- 
bation field ef when it is switched on. This makes the 
proposed internal mode mechanism for soliton ratchets 
quite general. 

Let us now briefly investigate the influence of a white 
noise on the kink ratchet dynamics. In Fig. |l2| we re- 
port the contour plot for the kink motion in the case of 
a bi- harmonic driver with m = 2. We see that the noise 
introduces disturbances on the profile but does not de- 
stroy the drift motion of the soliton. We checked that 
this property remains true also if we increase the am- 
plitude of the noise up to the kink-anti-kink nucleation 
limit. Moreover, the existence of the phenomena in pres- 
ence of noise shows the validity of the above mechanism 
also for the non-deterministic soliton ratchets. 



III. BOUNDARY EFFECTS ON SOLITON 
RATCHETS 

In this section we discuss the effects of the system 
boundaries on the ratchet dynamics. We have solved 
the problem for two types of boundary conditions: free 
ends Ux{0,t) — Ux{L,t) = and periodic boundary con- 
ditions u{0,t) = u{L,t) + 27m, n — ±1, where L is the 
length of the sample. The behavior of kink and anti-kink 
solutions does not differ for both cases except, of course, 
for the behavior at the boundaries. For free boundaries 
one can show by perturbation theory that the kink may 
be destroyed at the boundaries if its incoming kinetic 
energy is below a certain threshold. This can be easily 
understood since at the boundary the kink undergoes a 
large-amplitude oscillation which, in presence of damp- 
ing, expose the kink to larger dissipation (the oscillation 
will be damped and the kink will not be able to attain 
the —277 value and get reflected as anti-kink). The syn- 
chronization of the soliton motion with the external ac 
field, as well as, the fact that kink and anti-kink ratchet 
move is opposite directions, can allow sufficiently ener- 
getic soliton ratchets to overcome reflections in presence 



of damping. This is clearly shown in Fig. |l3| a where a 
kink-anti-kink ratchet reflection is shown. The possibility 
to overcome reflection, however, requires, to overcome a 
critical energy threshold which depends on the system 
parameters. In Fig. [Ts^ d we show the case in which 
the soliton ratchet is destroyed at the boundary. At low 
damping the collision of the kink with the boundary gen- 
erate oscillations which decay, after some time, leading to 
the destruction of the kink (increasing the damping the 
decay time quickly decreases). The possibility of over- 
coming reflection depends also on the relative phase of 
the driver and the internal mode oscillation on top of the 
kink (this dependence can be tested by shifting the initial 
incoming positions of the kink). In general, besides kink- 
anti-kink reflections and kink destructions, more compli- 
cated phenomena can arise. These including the possi- 
bility, for particular values of damping, amplitudes and 
driver frequencies, to sustain a standing breather at the 
boundary for very long time. Similar phenomena were 
also observed in Ref. | |2^ . 

To avoid the possibility of ratchet destruction at the 
boundaries one can recourse to periodic boundary condi- 
tions. In this case soliton ratchets, once established, will 
go on forever (we have checked this numerically for very 
long computation times). This opens the possibility of 
interesting physical applications as we will discuss at the 
end of the next section. 



IV. CONCLUSIONS 

In this paper we have considered a new way to produce 
a directed motion of a topological soliton in presence 
of damping, by using suitable ac drivers of zero mean. 
In particular we showed the possibility of establishing a 
ratchet dynamics for a kink of the damped sine-Gordon 
equation when a periodic force, consisting of two har- 
monic drivers, is applied. In contrast to previous works, 
the observed dc motion does not requires any asymmetry 
in the potential of system, this making the phenomena 
easily accessible to experimental situations. 

We also showed that a first order perturbation anal- 
ysis based on collective coordinates does not provide a 
complete description of the phenomenon. The reason for 
this failure was ascribed to the soliton-phonon interac- 
tion which, in a perturbation analysis, appears only at 
second order. In particular we showed that the soliton- 
phonon interaction manifests with the excitation of an 
internal mode on the soliton profile which, in presence of 
damping, can interact with its translation mode. This 
provides a basic mechanism to convert the energy of the 
ac field into direct motion for the soliton which is valid 
for a wide class of under-damped or moderately damped 
nonlinear systems. Numerical simulations of the sine- 
Gordon equation confirm the validity of the proposed 
mechanism. We also investigated the influence of bound- 
ary conditions on the ratchets dynamics in finite sine- 
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Gordon systems. For reflective boundaries we showed the 
pos sibility for sohton ratchet to overcome reflections in 
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tors: Noisy Transport Far from equilibrium. 



pres ance of dissipation (for periodic boundary conditions 



mat/0010237 (2000) 



the sohton ratchet, once estabhshed, wiU, obviously, go 
on forever). Finally, we showed that the phenomenon of 
soliton ratchet is robust enough to overcome the presence 
of the noise in the system. 

This results open the possibility of interesting appli- 
cations in different fields. In the context of Josephson 
junctions, for example, one can predict the existence of 
a zero field step (i.e. steps in the current-voltage char- 
acteristic related to the resonant fluxon motion in the 
junction) in absence of dc bias current and in presence of 
only bi-harmonic fields of zero average. This effect should 
be best observable in annular Josephson junctions where 
no boundary problems (kink destruction) exist. We also 
remark that "fluxon-ratchets" inducing zero field steps in 
Josephson junctions, have not yet been considered both 
theoretically and experimentally (we will investigate this 
problem in more details in a forthcoming paper). 

Similar transport phenomena can be predicted also in 
a variety of physical systems such as dislocations in crys- 
tals, spins waves in magnetic chains, etc. Adjusting the 
relative phase of the drivers one could control the direc- 
tion and the velocity of these excitations, this providing 
a way to control their dynamics. We hope that the re- 
sults of this paper will soon stimulate experimental work 
in these directions. 
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Figure captions 

FIG. 1. Profiles of tfie bi-harmonic driver for m — 2 (con- 
tinuous line) and m = 3 (dashed line) for parameter values 
El = 0.4, E2 = 0.26, uj = 0.25, 6 = 0. For a better compari- 
son, a time shift of 4.548 and of —2-k was respectively applied 
to m = 2 and m = 3 cases. 

FIG. 2. Time dependence of the velocity v{t) of the kink 
center of mass (thick dashed curve). The continuous line de- 
note the force profile, while the dotted-dashed and dashed 
curves represent the response to the single dark and light 
pulses, respectively. The parameters are the same as in Fig. 
hi with m = 2 and a = 0.15. 

FIG. 3. Trajectories of the kink center of mass derived from 
Eq. (jl^) for the cases m = 2 (continuous curve) and m = 3 
(dashed curve). The other parameters are: ui = 0.25, a = 0.1, 
El = 0.4, E2 = 0.26 and 6 = 1.61. 

FIG. 4. Average velocity (v) of the kink center of mass 
versus the the relative phase for the same parameters as in 
in Fig. ^. The continuous curve refers to the approximating 
function (t;> = 0.058 sin(6' - 0.8). 

FIG. 5. Contour plots of the velocity field Ux{x,t) with 
a = 0.12, Lo = 0.25, 6 = 1.61, in the cases {a)Ei = 0.4, 
E2 = 0.26 and (b)^! = 0.66, E2 = 0. 

FIG. 6. Contour plot of the velocity field Ux{x,t)of the 
anti-kink motion. System parameters are as for Fig. raa. 

FIG. 7. The averaged kink velocity as a function of damp- 
ing constant a with E2/E1 = 0.65 for different system param- 
eters. 

Panel (a): Ei = 0.4, 9 = 1.61; to = 0.11 (□), uj = 0.25 (o), 
LJ = 0.35 (•) and uj = 0.65 (o). 

Panel (b): to = 0.35, 9 = 1.61; Ei = 0.5 (*), Ei = 0.4 (•), 
El = 0.38 (□), El = 0.3 (o). El = 0.2 (+) and Ei = 0.1 (o). 
Panel (c): lo = 0.25, Ei = 0.4; 6 = (o), 9 = 0.8 (o), 6 = 1.2 
(□) and e = TT (*). 

The dashed lines show results obtained from numerical solu- 
tion of Eq. (0) for w = 0.11 and lj = 0.25 in (a); for Ei = 0.4 
in (b) and for ^ = in (c). 

The inset shows the dependence of the average velocity on 
the "delay angle" 9 for a = 0.15 and the rest of parameters 
as in Fig. |^. The dashed curve in the inset shows the fitting 
curve (see text for details). 



FIG. 8. Dependence of the averaged velocity (v) on damp- 
ing for parameters as in Fig. |^ but with different amplitudes: 
El = 0.4 (*), El = 0.45 (o) and Ei = 0.5 (o). The ratio 
E2/E1 has been kept constant, E2/E1 = 0.65. 
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FIG. 9. The averaged kink velocity as a function of the 
driver amphtude Ei with E2/E1 = 0.65, d = 1.61. The rest 
of the parameters were fixed as: 

ij = 0.11, Q = 0.15(*); io = 0.35, a = 0.15 (o); lj = 0.35, 
a = 0.1 (o); = 0.35, Q = 0.05 (□). The inset shows the 
log-log dependence. 

FIG. 10. Contour plot of the velocity field Ux{x,t) for a 
sine-Gordon kink driven by a bi-harmonic force with m = 3. 
System parameters are as for Fig. ^ except a = 0.15 and 
e = TT. 

FIG. 11. Kink profile and dynamics for a — 0.15, uj — 0.35, 
El = 0.5, E2 = 0.325, e = 1.61 (a,b) and = 1.61-n = -1.53 
(c). The profile has been computed at the time moment 
t — 600. The dashed line shows the profile of the kink af- 
ter one period of the external drive. 

FIG. 12. Countour plot for the displacement field u{x,t) 
for parameters as in Fig. la and noise amplitude D = 0.1. 

FIG. 13. Contour plots of kink (a) reflection and (b) anni- 
hilation at the boundary, the system parameters are uj — 0.11, 
El — 0.4, E2 = 0.26, 6 = 1.61. The damping parameter was 
a = 0.08 for the case (a) and a = 0.12 for the case (b). 
The inset shows the anti-kink profile after reflection at time 
t = 3200. 
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